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I. INTRODUCTION 



According to Ronveaux (1995) the Heun equation (HE) is the most general second order 
linear ODE of the form 

d 2 y /l + 2«i l + 2a 2 . 1 + 2a 3 \ dy af3z-q A _ /lN 

where a G C\{0, 1}, ai, a 2 , a 3 , a, /3, q A are complex arbitrary parameters, and 0, 1, a, and 
oo are regular singularities with exponents {0, — 2ai}, {0, — 2a 2 }, {0,-20:3}, and {«,/?}, 
respectively. Equation ([1]) has been originally constructed by the German mathematician 
Karl Heun (1889) as a generalization of the hypergeometric equation (HYE). In order to see 
how HE degenerates to the HYE we can first multiply ([I]) by z(z — l)(z — a), then we set 
a = 1, and q A = a/3, and finally we take out a factor (z — 1), leaving the HYE in its standard 
form. Hence, we can always think to a HYE as a degenerated equation of Heun's type. To 
underline the importance of ([1]) we recall that it contains the generalized spheroidal equation 
(GSWE), the Coulomb spheroidal equation, Lame, Mathieau, and Ince equations as special 
cases. The fields of applications of the HE in physics are so large that it is not possible 
to describe them here in detail. However, a review of many general situations relevant to 
physics, chemistry, and engineering where the HE occurs can be found in Ronveaux (1995) 
(pp 341). Here, we will show which role plays the HE in quantum field theory in curved 
spacetimes. 

To understand the motivation underlying the present work we give a short review on studies 
concerning exact solutions of the Teukolsky equation in some black hole geometries. In 
the 70's, and 80's we find a large number of publications regarding the angular equation 
obtained after separation of variables from the Teukolsky wave equation on Kerr manifolds. 
See, for instance, Press, and Teukolsky (1973), Breuer et al. (1977), Fackerell, and Crossman 
(1977), Leahy, and Unruh (1979), Chakrabarti (1984), and Seidel (1989). According to these 
references we will also name the solution of the angular equation as spin- weighted spheroidal 
function (SWSF). The radial equation has been investigated by Bardeen, and Press (1973), 
Page (1973), Lee (1976), Arenstorf et al. (1978). The common picture emerging from 
all previous studies is that the radial equation cannot be in general related to any known 
differential equation of mathematical physics. This view changed with the work of Blandin 
et al. (1983). They showed that the SWSF's may be obtained by means of an elementary 
transformation from Heun confluent functions. Three years later Leaver (1983) proved that 
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the radial, and angular parts of the Teukolsky master equation (TME) in the Kerr geometry 
are generalized spheroidal wave equations. Finally, Suzuki et al. (1998) showed that the 
radial, and angular part equations arising from the TME in the Kerr-Newman-deSitter 
metric (KNdS) after separation of variables are Heun equations. It is interesting to observe 
that if we let the cosmological constant go to zero (i.e. the KNdS geometry goes over to 
the Kerr-Newman metric) their HEs become a confluent HE which, in turn, coincides with 
the GSWE given by Leaver in 1983. Hence, the following question arises quite naturally, 
namely: is it possible to reduce the TME in any physical relevant type D metric to a HE? 
To conclude this short review we cite what Wu, and Cai wrote in 2003: "it is not clear until 
now whether the generalized Teukolsky equation in the general type D vacuum backgrounds 
with cosmological constant can be transformed into a Heun equation." 
Our paper is organized as follows: in Sec. [Til we shortly present some results due to Kamran, 
and McLenaghan (1987) concerning the separation of the TME in any type-D background. 
In Sec. lITTllVHl we show that the TME can be transformed in any physical relevant type D 
metric into a HE. 



II. BACKGROUND 



Let denote the class of algebraically special Petrov type D vacuum metrics with 
cosmological constant. According to Thm. 2.1 in Kamran, and McLenaghan (1987) there 
exists a system of local coordinates (u,v,w,x) in which such metrics can be written as 



ds 2 = 2 {e l e 2 - 9 3 9 A ) 

with a symmetric null tetrad {9 , 6 2 , 8 3 , 8 4 ) given by 



o 1 
e 2 



V2T(w,x) 

y f Z(w,x) 

V2T(w,x) 

y f Z(w,x) 

V2T(w,x) 



fW(w) 
Z(w, x) 

W{w) 
Z(w, x) 

X(x) 
Z(w, x) 



(ei du + m(x) dv) + 
(ei du + m(x) dv) — 



dw 



(e2 du + p(w) dv) + i 



g 2 W(w) 

fdw 
g 2 W(w) 

dx 



X(x) 



Z(w,x) := eip(w) — €2m(x), g :- 



(2) 



where all functions are real-valued and ei, 62, and / are constants such that e\ + e 2 , 7^ 0. 
Depending on whether fW 2 (w) is positive, negative or zero the metric (j2J) possesses a two- 



parameter abelian group of isometries whose orbits are timelike, spacelike or null at a given 
point, respectively. By integration of the Einstein-Maxwell field equation it results that the 
general solution srf* in the class @ Q can be specified as follows (Thm. 2.2 ibid.) 



m{x) 

p(w) 
T(w, x) 
fW\w) 
X 2 (x) 

K 2 



b 2 cos 7, e 2 = sin 7, 
cV + b 2 k 2 + 

b 2 c 2 w 2 + f + k2 ( b2 - e ^ 



;i-6 2 cos 2 7 ) 



e 2 — 2c£x, 



cos 2 7 



cos 7 + 2b 2 ckw, 



a(cw cos 7 + k)(cx sin 7 + £) + 1, 
c 2 6 4 5f 4 u; 4 cos 2 7 + C/3W 3 cos 7 + f 2 w 2 + fiw + f , 
c 2 g4X 4 sin 2 7 + c/tix 3 sin 7 + /t 2 x 2 + + g , 
acfi cos 7 — 2akf 2 + 3ak 2 f 3 + A{£ — ab 4 k 3 )g4, 

3ac£/i cos 7 - (1 + 6ak£)f 2 + Sk(l + 3aki)f 3 + 6[£ 2 - 6 4 A; 2 (1 + 2ak£)}g A 



where fo, fi, f 2 , fz, go, <?i, 94, a , b, c, k, £, and 7 are real parameters satisfying the relations 



acgi sin 7 - 3a 2 ctfi COS7 + 2a£(l + 3ak£)f 2 - [1 + 3ak£{2 + 3ak£)]f 3 



+ A[b 4 k - at + 3ab 4 k 2 £{l + ak£)]g A = 0, (3) 



c 2 (g sin 2 7 - b 4 f cos 2 7) + c[(2a£ 3 + b A k)h cos 7 - lg x sin 7] - (& 4 fc 2 + £ 2 + Aak£ 3 )f 2 

+ (b 4 k 3 + 3k£ 2 + 6ak 2 £ 3 )f 3 + (3£ 4 - b s k 4 - 6b 4 k 2 £ 2 - 8ab 4 k 3 £ 3 )g 4 = 0, (4) 

and are restricted to a range such that (j2J) is non-singular with signature minus two. More- 
over, the cosmological constant A is expressed in terms of the above parameters by 

A = -3[a 2 c 2 / cos 2 7 - a 2 ckf x cos 7 + a 2 A; 2 / 2 - a 3 k 3 f 3 + (1 + a 2 b 4 k 4 )g 4 }. (5) 

Let srf* denote the subclass of solutions in srf* obtained by setting b = 1, c = v2, k = £ = 
and 7 = 7r/4. If in addition a — f — 1, such solutions recover the vacuum case with 
cosmological constant of the seven-parameter family of Plebanski and Demianski (1976) 
containing the Kerr-Newman-de Sitter metric as a special case. 

Let be the subclass of solutions in stf* such that a = 0, b — c — 1, £ = and 7 = 7r/2. 
If in addition / = 1, SS Q _ reduces to the vacuum case with cosmological constant of Carter's 
B_ (1968). By we will denote the subclass of solutions in srf* obtained by setting a = 0, 
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b = c = 1, and k = 7 = 0. If we let / = 1, such a solution becomes the vacuum case with 
cosmological constant of Carter's B + . The Carter's B± describe all non-accelerating type 
D metrics in a coordinate system in which the components of the metric and the Maxwell 
field are rational functions. 

Let ^* denote the subclass of solutions in sd* obtained by setting a = 1, b = 0, c = \/2, 
k = i = and 7 = ir/A. If in addition / = 1, c ta* reduces to the accelerating ^-metric of 
Levi-Civita (1918). 

Finally, let ^ 00 denote the subclass of solutions in &/* obtained by setting a = 0, b = 1, 
c = 0, k = 0, £ = 1 and 7 = 0. If in addition / = 1, becomes the Robinson-Bertotti 
solution (1959). 

Following Kamran, and McLenaghan (1987) the Teukolsky equation can be written in a 
compact form by introducing a spin parameter s which can assume the values ±2, ±1, and 
±1/2. For s = 2, 1 and 1/2 we have 

[(£> - (2s - l)e + e - 2sp - p)(A - 2sj + /x) - {5 + W - a - (2s - l)/3 - 2st) [5 + tt- 2sa) 

-2(a-i)( S -l)* 2 ]*, = 0, (6) 

and for s = —2, —1 and —1/2 

[(A - (2s + 1)7 - 7 - 2s/i + p) (D - 2se - p) - (5 - r + - (2s + l)a - 2six) (5 - r - 2s/5) 

-2(s + i)(s + l)^ 2 ]$ s = (7) 

where ^2 is the non-zero Newman-Penrose component of the Weyl tensor and the <3> s 's are 
defined in terms of the field components as given in Table 1 (pp 286 ibid.). According to 
Thm. 3.1 (ibid.) for all solutions in the class £>o, and for all s = ±2, ±l,±l/2 equations 
(EJ), and (j7|) possess a separable solution of the form 

where r, q are arbitrary real constants, dB = Z^ 1 (eim dw + e 2 p'dx), and Q s (w,x) = 
G s (w)H s (x). Moreover, (jSJ), and © separate into the pair of decoupled ODE's 

D ws L ws G s (w) + f s {w)G s (w) = \ s G s (w), (8) 
D xs L xs H s (x) - g s (x)H s (x) = \ s H s (x) (9) 
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where X s is a separation constant, the functions f s , and g s are given in Table 1 (pp 290-291 
ibid.), and 



D u 

L u 
L, 



W— + 



f(l + e(s))-j(l-e(s)) 



W 



dw 1 + f 2 

. v d . , ,t\q - mr ... dX 

-iX — + ie{s)-?— + i( s - 1 — 

ax X dx 



pr-e 2 q , . ,(W 
+ (1 - |s|)- 



1 + / 2 



[(1 + cW) - f (1 - e(s))} 



W 



W 



dW 

dw ' 



d e\q — mr dX 

lX Tx +ie{s) ^^ + M ^ 



where e(s) is a sign function such that e(s) = +1 for s > 0, and e(s) = —1 for s < 0. 



III. THE METRIC 



b = 1, c = V2, Jfe = £ = , 7 = tt/4 



In this case we have 



72 



ei = e 2 = — , m(ae) 



— V2x 2 , p(w) = V2w 2 , T(w, x) — awx+ 1, 

4 

n=0 i=l 

4 

X 2 (x) = 5f 4 a; 4 + a/ix 3 - / 2 a: 2 + #i£ + g = g A - Xi). 

i=l 

Moreover, for i,j — 1, • • • ,4 let w« and Xj denote the i-th root of the polynomial equations 
fW 2 (w) = 0, and X 2 (x) = 0, respectively. Throughout this section we shall assume that 
Wi 7^ Wj, and Xi ^ Xj for every i ^ j. Furthermore, in the present case (jSD reads 



O 2 /0 + 77 ) • 



#4 = • 

Making use of the expression for \I/ 2 given by (2.7g) in Kamran, and McLenaghan (1984), 
the functions f s (w), and g s (x) entering, respectively in (jSj), and Qj are computed to be 

f.(w) = -(1 - 2|s|) [2^4(1 - \s\)w 2 + (2iV2re(s) + ag x {\ - |s|)) w 
<? s (a;) = -(1 - 2|s|) [2^4(1 - \s\)x 2 - (2V2re(s) - af^l - |s|)) x 

Hence, the equation for G s (w) becomes 

s(/^§)+r(«)§ + 0.(«)G. = o 



(10) 



with 



f 2 - 1 

F(w) := iV2 J — -{2rw 2 - q), 



Q a {w) := \ s +2g 4 (l-\s\)(l-2\s\)w 2 + 



P + l 



2i\[2r ( -2s + 



J 2 -l 
P + l 



+ a^(l-H)(l-2|s|) 



w+ 



f\ s \ f WW " + (1 _ | S |)(^') 2 ) + zV2 5 (2to 2 - g)^ + 21 (W ~ g) ' 



W (l + P) 2 W 2 
Equation ( jTUl) can be further simplified. To this aim let us make the substitution 



G s (w) = e^V»- 



If we require that h' 



-T/(2fW 2 ) we obtain the following equation for (p s (w) 

+ R s {w) <p s = 



d_ 

dw 



dw 



(11) 



where 



R s {w) := \ s + 2^(1 - |s|)(l - 2|s|)«; 2 - C s w+ 

f\s\ (WW" + (1 - |,|)(^') 2 ) + i V2s(2rw 2 -q)^ + 

with 

C s := 4iV2sr - ag^l - |s|)(l - 2|s|). 
Let us introduce the following functions 

a±q(w) := 2rw ± q, f s (w) := 2g 4 (l - |s|)(l - 2\s\)w 2 - C s w + X s , 



and let us define constants 



c i 1 := 94 ~ Wj), i = 1,2,3. 



3=1 

By means of the homographic substitution 

w — W\ W2 — W4 

z = 

W — W4 W2 — Wi 

mapping the points wi,w 2 ,W3,W4, 00 to 0, 1, z$, 00, Zoo with 



w 2 — W4 
w 2 - Wi 



zs 



W3 — W4 



(12) 



(13) 



equation (TTTT) becomes 



d 2( Ps sd<p s ~ , . 



(14) 



with 



/ x 1 1 1 
P( Z ) = - + + 



Z Z — 1 Z — Zs Z — Zoo 

B x B 2 B 3 2 A x A 2 A 3 A- 



+ 



z 2 (z - 1) 2 (z - Z S ) 2 (z - 2oo) 2 Z Z — l Z- Z S Z- Z a 
where for i = 1, 2, 3 

3 



A r 



B, 



1 



^oo(w 4 - Wi) 



(2\s\ 2 - 3\s\ + 2) J2 w i + C 1 - I 5 !)! 1 - 2|s|)^ - (2 - |s|)(l + 2|s|K 



i=l 



S _ . CiCT-q^Wi 



2\ \ 2 



2 \/2(u>j - ^4) 



w 2 - Wi 



Ai = Qs(wi), A 2 = c 2 - 

Zoo W4 - W\ 



Q s (w 2 ), A 3 = c 3 



w 3 - Wi 
z s {w A - Wi) 



9s(w 3 ), 



and 



Q s {wi) = ^\s\p(wi) - f s {wi) + i V ^ S cr+g (wi(w 4 - 2wi)) + 
2 Wi — W4 



Wi — W4 



a -q( w i) 



+q 



3=1 



n 



W 1 



\ 



p(wi) :=2(2|s|-3)^ 2 + (4-3|s|)^^^ + (|s|-2)w 4 ( ^w j -2w i +2(|s|-l) JJ 

/ 



3 = 1 



W j 



3=1 

3^i 



(15) 



If we make the F-homotopic transformation 



<p s {z) = z a \z - \T\z - z s ) az {z - Zoo )$ a {z) 



and require that of = —Bi, then (fl4"|) becomes 



-— + P s (z)— + Q s \z)y 8 = 

dz z dz 



(16) 



with 



1 + 2«i , 1 + 2ct 2 , 1 + 2a 3 - Ai A 2 A 3 A^ 



+ 



+ — — Qs(z) = — 



z—l z — Zs 



+ — — + 



z z — l z — Zs z — z. 



s 



Now, a direct computation gives 



(l-N)(l-2|s|) 

Zoo(w 4 ~ U>l) 



Wi + — 

94 



J=l 



:i-|a|)(l-2| a |) 

^oo(w4 - Wi) 



(17) 
(18) 



i=l v '\ i=l 

Taking into account that for the metric srf* equation ([3]) reduces to agi — f 3 — 0, and 
replacing /3 by agi in the equation /W 2 = it can be checked that the following relation 
holds 

4 

agi 



Hence, the coefficients (iTTj) . and (1T8"]) are zero, and (|T6l) reduces to an Heun equation with 
terms aft and accessory parameter q A (according to the notation in ([1])) given by 



a(3 = -{z 3 + l)Ai - z s A 2 - A 3 = z s A 3 + A 2 , q A = -z s A 1 



(19) 



where 



Ai = Ai — an — a 2 (l + 2ai) 

-r / a 2 + a 3 (l + 2a 2 ) 

A 2 = A 2 + ai + a 2 (l + 2ai) - — ^- - 

zs — 1 

a 2 + a 3 (l + 2a 2 ) ai + a 3 (l + 2a x ) 



1 



A 3 = A 3 + 



+ 



*S - 1 ^5 

Concerning ([9]) we find the following equation for H s (x), namely 



z — 1 



Z'Y- Za 



(20) 
(21) 
(22) 



d_ f x2 dH s 
dx V c^c 



-A s + 2^(1 - |s|)(l - 2|s|)x 2 + C s x+ 
(XX" - (|s|- 1)X' 2 ) - ^2 S (2 ra 2 + - + 



# s = (23) 



with 



C s := 4^2rs + (1 - |s|)(l - 2|s|)a/i. 
Let us introduce the following notation 

U{x) = 2g 4 (l - \s\)(l - 2\s\)x 2 + C s x - A 6 
3 

Ci ~ l = g4 Y[(%i -Xj), i = 1, 2, 3. 
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By means of the homographic substitution 



X — X\ x 2 — x 4 
z = 

X — £4 X2 — X\ 

mapping the points X\, X2, x 3 , X4, 00 to 0, 1, xs, 00, Xoo with 



(24) 



_ x 2 - Xi _ X 3 - Xi 

Xqq , Xs %oo 



x 2 - X\ 



x 3 - x 4 



our equation (|23|) becomes 



d 2 H, 



(25) 



(26) 



with 



*P.(*) = - + ~^-r + ' 



Z Z — 1 Z — Xs Z — Xqo 

B x B 2 B 3 2 

" + 1 T7T + 



Ai A 2 A 3 Aoo 
+ — + — — + — + 



Z 2 (z — l) 2 (z — Xs) 2 (z — Xoo) 2 z z — 1 z — Xs z — x 



where for i = 1,2,3 



U = — ^ r I (2\s\ 2 - 3|a| + 2) y> 4 + (1 - |a|)(l - 2\s\)^ - (2 - | a |)(l + 2\s\)x 4 ) , 

X00KX4 - xi) \ j~f g A I 



i=i 



Bi = - 



and 



£ + Cj a +q (x 2 ) 
2 y/2(xi - x A ) 



A x = — — Qs{xi), A 2 = c 2 X2 _ Xl Q s {x 2 ), A 3 = c 3 Xs Xl Q s (x 3 ), 

Xqq X^. X\ Xq\X^ X\) 



Qs(xi) = ^\s\p(xi) - f s (xi) y^s_ ^ ( x .( Xa _ 2xi)) 

Zi Xi X4 



X i X 4. 



-a+gix 2 



<J- q {x 2 ) - 2Xi<J_ 



3=1 



n 



p{x t ) := 2(2|s|-3)x 2 + (4-3|s|)x i ^x j + (|s|-2)x 4 

3=1 

If we make the F-homotopic transformation 



/ 3 \ 
^Xj-2xi +2{\s\-l)J\x j . (27) 

/ 



3=1 



# s (z) = z 3l (^ - 1) S2 (^ - ^) S3 (^ - x^Hsiz) 
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together with the requirement that a 2 = —B^, then fl56|) becomes 

d 2 H s - . ,d# s 



+ <p s (z)-^ + £2 s (z)if s = (28) 



with 

^ s (-2) = 1 r" H ' &s{z) = 1 r H h 



oo 



Z Z — 1 2 — " ' ' Z 2! — 1 Z — Xs Z — X c 

Now, a direct computation gives 

;i-H)(i-2H) a/, 



a.= ^"7 A1 ~7M E" + — ■ < 29 > 



^ (1-|.|)(1-2|.|)/A gA\ 

Going back to the polynomial equation X 2 (x) = it can be verified that 

4 t 

Eaji 

1=1 y4 

This implies that (|29|) . and (|30|) are zero, and (128]) becomes an Heun equation with 

a(3 = x s 2l 3 + 2l 2 , g A = -x s 2li (31) 

where the can be directly obtained from the formulae ([2"U1) - (T2"2"1) for the coefficients A{ 
by means of the formal substitutions Ai — ► 2lj, Aj, — > A^ ot{ — ► 5?j, z$ — ► Xs and 

IV. THE METRIC ^° : b = c = 1, a = £ = , 7 = vr/2 

In the present case we have 

ei = 0, e 2 = 1, m(x) = — (x 2 + k 2 ), p(w) = 2kw, T(w,x) = 1, 
fW 2 (w) = f 2 w 2 + fxW + f = f 2 (w - wx)(w - w 2 ), 

4 

X 2 (x) = g A x 4 + k 2 x 2 + g x x + g = g 4 JJ(a; - x { ), k 2 = -f 2 + 3kf 3 - 6k 2 g 4 . 

i=l 

Moreover, for i,j — 1, • • • , 4 let w±, w 2 , and Xi denote the roots of the polynomial equations 
fW 2 {w) = 0, and X 2 (x) = 0, respectively. Throughout this section we shall assume that 
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Wi 7^ w 2, and Xi ^ Xj for every i ^ j. Finally, let us recall that in the present situation 
flU) and (jSJ) read, respectively 

4*34-/3 = 0, 9o~k 2 (k 2 g,-kf 3 + f 2 ) = 0, A = -3g 4 . (32) 

Since the term x 3 is not present in the expression for X 2 it can be easily checked that the 
roots of the polynomial equation X 2 (x) = satisfy the condition x\ + x 2 + x% + £4 = 0. 
Taking into account the expression for \&2 given by (2.8g) in Kamran, and McLenaghan 
(1984), the functions f s (w) and g s (x) entering, respectively in (El), and (jHJ) are computed to 
be 



f s (w) = 0, 

. IW , . j ... . ,. \ (X 2 (x))' QX 2 (x) + C A (k + ix) 
g.{x) = -2ire(s)(l - 2\s\)(k + ix) + (1 - | a |)(l - 2|«|) ^±^~ 3 ( fc _ tg) 2 " 

with C*a := —2k(3f 2 + 4AA; 2 ) + 3igi. The equation for G s (w) becomes 



with 



f -I 

T(w) := 2i——(2rkw - q), 
/ + 1 



Q s (w) := \ s - (1 + e{s) - f{l - e{s))) + f\s\ (ww" + (1 - \s\)W' 2 ) + 



2is(2rkw - a)— + 4/ ^ rkw - & 
2zs(2rkw D w + {l + f2)2 w2 ■ 

In order to simplify (1331) we make the transformation 

G s (w) = e^V-H 

requiring that ft,' = — Y / (2fW 2 ). Hence, we end up with the following equation for <p s (w) 

d ( j. TJ/2 dip s 



^{ fW ^) +R -^- = ° < 34) 

where 

,2\ W' (2rkw - q) 2 



/ 11 i2\ W 

R s (w) := X s - 2irke(s) + f\s\ [WW +(l-\s\)W J + 2is{2rkw - q)— + 
By means of the transformation 



f W * 



z = (35) 
w 2 - Wi 
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mapping the points wi, w 2 , 00 to 0, 1, 00 our equation (I34I) becomes 



with 
where 



- wr + P.(z)-^ + Q.iz)<p. = 



p f \ 1 1 n n Sl ^2 A 2 
P a U) = -H r , = — + 7 — H h 



2; 2! — 1 



s . 1rkw\ — q 

— ~t~ 1 

2 h(w 2 -wi) 



Bo 



z 2 (z — l) 2 2; 2 — 1 



s 2rkw2 — q 
i- 



2 f2(w 2 -w 1 ) 



h 2 ' / 2 > 2 - Wl )2ll 

If we make the F-homotopic transformation 

ip s (z) = z a ^z-ir^ s (z) 

with the requirement that a 2 = —Bi then (1361) reads 

(z) + P»£C(z) + Qs{z)y s {z) = 



(36) 



A s + 2irk(e(s) - s) |s|(|s|-2) 2 A/ , x 



-M- 



(37) 



with 



3 M _ l + 2ai l + 2a 2 - _ (Ai + A 2 ) ^ + a 2 + «i(l + 2« 2 ) - A x 



z z — 1 - - > • ^ _ 

Since A 2 = — A\ it follows that (1371) is an HYE. Notice that an HYE can be always thought 
as a degenerate Heun equation since it can be recovered from (JT]) by setting for instance 
Zs = 1 and = or z s = q = or 1 + 2a 3 = and q = z$ ct[3. Concerning the ODE (J9]) 
for HJx) we find 



d 

dx 



( x2 ^) + [~ (As ~ ^^X 1 " 2 I S D) + 4rsx + I s ' { xx " + ( x _ l s D x ' 2 ) - 



2r,(a; 2 + fc 2 )^-r 2 (3:2 +f )2 - 



(1-N)(1-2H) 



Let us introduce the following function 



(X 2 )' _ 6X 2 + iC A (x-ik) 
x + ik 3(x + ik) 2 



H s = 0. (38) 



a±(x) := x ± k 2 
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Moreover, let q be defined as in the previous section. By means of the homographic substi- 
tution fl24l) equation (1381) becomes 



with 



d 2 H* . . dH* . . . „ 

_£ + ^ (2) _l + Q , (2)ft = 



¥.(*) = - + 1 



z z — 1 2; — X5 2 — x c 



(39) 



jQ 5 (z) = -^ + - — + - — + - - + -, ^TT-I- 



(Z - l) 2 (Z - X S ) 2 (Z - Xoo) 2 (Z - X k y 

Ai A 2 A 3 joo Ife 



2; 2 — 1 -2 — X5 Z — Xqo Z — Xk 



5 >£fc • -^c 



where for i — 1, 2, 3 
1 



^Oo(^4 ^l) 



^^^^^ ^2 3fcCi 



4 , 



.1=1 



s re, a^(x 2 ) ^ 2 
7; H _ I > ^ 



'x 4 + ik 



kC A 



-2(i-| S |)(i-2| S |)(i+^-An^+^) -1 )' 



J'=l 



7 ci c ~ x 2 - x x ~ c 3 (x 3 - Xi) 

= Ba{xi), A 2 = c 2 Q,{z 2 ), A 3 = — s (x 3 ), 

Xoo £4 - Xi X 5 (X 4 - Xi) 



(1- | g |)(l-2|s|)(x 4 + ife) 
(xi - x 4 ) 



.C 



(k) + i^p{k) \[{x j + zfc)- 1 ) l[(x t + ik)-\ 



i=i 



i=l 



with 



#4, 



2rs 



2^2 



Qs(xi) = X s - 2irke(s)(l - 2\s\) + ^\s\p(xi) + 



cr + (XjX 4 J o + l^ 



(T_(x 2 ) Xj — 2XjO"_ 



n 



(1 - |s|)(l - 2|s|) / x,j - x 4 _ i ^ Xj -ik \ 
Xj + ik V Cj 3 Xi + ik J : 



r(k) := ^3x 4 — ^^^j^ ^ 2 + 2i(xiX2 



+ X1X3 — X2X 4 — X3X4 — X]X 4 + X2X 3 )k+ 

3X1X2X3 — X2X3X 4 — XiX2X 4 — XiX3X 4 , 



p(k) := 3k 4 — I Xj + 5x 4 J k 3 + ( XiX 2 + x 2 x 3 + X1X3 — 3x 4 X; L J k 2 + 



i=i 



i(xiX2X 4 + XiX 3 X 4 + X2X3X 4 — 3X1X2X3)/^ — XiX2X3X 4 , 
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and p(xi) given by (|27|) . Inserting the roots of the polynomial equation X 2 (x) = into the 
expressions for Bk, and Ak we checked with the help of Maple 9.5 that Bk = Ak = 0. Let 
us work with the F-homotopic transformation 

H s (z) = z^(z - lf^z - x s f 3 (z - Xoo )H a (z). 

Moreover, we require that af = —B^. Then (1391) becomes 

H" s {z) + $ s (z)H' s (z) + &s(z)H s (z) = 

with 

ft / x 1 + 25, 1 + 252 1 + 253 a , , 2li , 2l 2 B 3 
<p s (z) = 1 — H , U 8 (z) = 1 + 



(40) 



z z — 1 z — X s ' " ' £ £ — 1 ' z — x_ 

Taking into account that the 2lj's can be directly obtained from the formulae (|20[) - 

the coefficients Aj by means of the formal substitutions Aj, — > 2lj, Aj — > Aj, 

£s — > xs and — > Xoo, and that Ak = implies that 



for 



we obtain 



1 1 1 
+ + 



i=l 



i=l 



(1_ | s |)(l_ 2 | s |)(x 1 -x 2 ) 



(x A - x 2 )(x 1 - X A ) 



3x 4 - 22 Xi + 



i=l 



TT (xi + iky 1 T(k) - i—p(k)(x 4 + zA;) TT(^ + ^r 1 



(l-H)(l-2|s|)(x 1 -x 2 ) 



M(k) 



with 



(x 4 - x 2 )(xi - x 4 ) 

3 3 

Af(fc) := 3^ ~ E Xi + ( T ( fc ) + ( X4 + ^M*)) 11^ + 



(41) 



8=1 



T(/c) := | ^^Xj(2x4 — Xj) — 3^4 J k 2 + 



i=i 



( 



\ 



Ax 4 



z>+n 



fc + X 



I ,=1 

\ Ml 



1 = 1 

Ml 



/J 





i- 


=1 






3 








2^+E^ 2 








> 




/ 3 


3 ^ 




3 






+n^ 


\ Mi 


i=l 1 

Mi / 







.1=1 
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Inserting the roots of X 2 = into (|4"T1) it can be verified that M(k) = for all real k. Hence, 
we have reduced PU|) to an Heun equation with terms a/3, and qA given by fl3Tl) . 



V. THE METRIC <M\ : b = c=l, a = k = j = 



In the present case we have 



€i = 1, e 2 = 0, m(x) = — 2£x, p(w)=w 2 + £ 2 , T(w,x) = l, 



fW 2 (w) = g,w 4 + f 2 w 2 + f lW + f = g, l[(w - 



i=l 



X 2 (x) = k 2 x 2 + gix + g = k 2 (x - xi)(x - x 2 ), k 2 = -(f 2 + 2£ 2 A) 

where for % = 1, • • • , 4 Wi, and x±, x 2 denote the roots of the polynomial equations fW 2 (w) = 
0, and X 2 (x) = 0, respectively. Throughout this section we shall assume that Wi ^ Wj for 
every i ^ j, and X\ ^ x 2 . Finally, since the term w 3 is not present in the expression for fW 2 
it can be easily verified that the roots of the equation fW 2 (w) = satisfy the condition 
W1+W2+W3+W4 = 0. Taking into account the expression for ty 2 given by (2.9g) in Kamran, 
and McLenaghan (1984), the functions f s (w) and g s (x) entering, respectively in ([8]), and ([9]) 
are computed to be 

" fW 2 ' QfW 2 -C A (w + i£ 



f s (w) = -2ire(s)(l - 2\s\)(w + it) - (1 - |s|)(l - 2|s|) 
9s{x) = 



w 



3(w — it 



with C A := -2i£(3g 2 + AM 2 ) + 3/i. Concerning the ODE © for G s (w) we find 

aw V aw } aw 



(42) 



with 



f 2 - 1 

T e (w) = 2ir— — -(w 2 + £ 2 ) 



f 2 + l 



Q s ( w ) = X s - 2e(s)(l - 2\s\)r£ + 2ire(s)(l - 2|s|)w - 



2zr 



1 + f 2 
IT' 4/r 2 



[1 + e ( a ) _/»(!_ e(a ))] w+ 

(Viy" + (i - |s|)iy' 2 ) + 2tr S (w 2 + e 2 )^r+ a " ^ 7 + 

(i-H)(i-2H) 



+ /2)2 V^ 2 

/W 2 ' 6/H/ 2 -C A (w + ^ 



w — i£ 



3(w — U 
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By means of the transformation G s (w) = e h ( w ^ip s (w) with the requirement that h 
-V e /(2fW 2 ) equation fl42]) simplifies to 

d f j, rxr2 d(p s 



^l fW ^i +R -^- = ° < 43) 



where 

~ / 9\ W r' 2 (w 2 4- f 2 ) 2 

R s ( w ) = A s - 4zr SW + f\s\ [WW + (1 - ) + 2irs(w 2 + f)— + 1 fTT T 9 } 4 

V / W jW 2 

(i-M)(i-2M) 



/W 2 ' 6,/W 2 - C A (u; + U 



w — il 3(w — il) 2 

with 

\ s := A s - 2e(s)(l-2\s\)r£. 

For notational purposes let us introduce the function a±(w) := w ± £ 2 . By means of the 
homographic substitution ( fl2|) equation ( |43l) becomes 

^ + ^.(^+5.(^ = (44) 

with 

/ x 1 1 1 2 



Q ( ) ^ 1 + ^ 2 + ^ 3 + ^ + ^ + 

Z 2 (z - l) 2 (z - Z S ) 2 (Z - Z^) 2 - Zi) 2 

A 1 A 2 A 3 A t Wl 
1 rH 1 1 , z*:=2oo 

Z Z — 1 2 — Zs Z — Zoo Z — Zi W4 

where for % = 1, 2, 3 



Aoo = T— \ I ^ W i ~ 3W 4 



,1=1 
,2\\ 2 



A Cl / \ A W 2 -Wi C 3 (w 3 - Wl) 

Al = fls(Wl), v4 2 = C 2 0sW, ^3 = 7 r Q S {W 3 ), 

Zoo Wi — Wx Z S {W 4 -W 1 ) 



Wj — n 



17 



with 



g s (u7i) = -A. s + - spiflj a + {wiW 4 

2 Wi — W4 



2^,2 



2g 4 cfr 



-(T+{wj)- 



( , \ 



o- w; 



! ) Wj — 2wi<j- 



n 



W ; 



W{ — U>4 

(l-|s|)(l-2|s|) f Wi -w 4 



w, 



i pi Wi + i£ 

o^A ~ 



Wi 



?(£) : = ^3w 4 - y^WjJ 



?>WiW2Wz — W2W3W4 — W1W2W4 — W1W3W4, 



WIW 2 



+ w 2 w 3 + W1W3 - 3u> 4 



Wi 



.1=1 



i=l 



i(WiW 2 W 4: + W1W3W4 + W2W3W4 — 3WiW2W 3 )£ — WiW2W 3 W4, 

constants defined as in Section [TTTj, and p(wi) given by (|T5|) . Inserting the roots of the 
polynomial equation /W 2 (tu) = into the expressions for B^, and At we checked with the 
help of Maple 9.5 that Bf — Af — 0. In order to reduce (jUJ) to a Heun equation we just 
need to introduce the F-homotopic transformation 



<p a (z) = z a \z - \) a \z - z s T{z - Zoo )lp s {z). 

with exponents ai such that a? = —Bi, and to proceed as we did for the equation (140]) in 
Section HVl Concerning the equation for H s (x) we find 



A. (x 2 — 
dx V dx 



+ 



X 

-X s + 2e(s)£r - 2s(2£rx + q)—+ 

X 



\s\ (XX" + (1 - | S |)X' 2 ) 



,2\ {2£rx + qf 



X 2 



H s = 0. (45) 



By means of the transformation z = (x — X\)/{x2 — X\) mapping the points Xi,X2,oo to 
0, 1, 00 our equation (1451) becomes 



^+?,(z)^ + fl,(z)ff, = 



(46) 



with 



1 1 

+ 



2; z — 1 



aw 



-^1 j #2 



+ — + - 



2; — 1 
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where 



. s 2r£xi + q \ 2 ~ / s 2r£x 2 + g x 

&\ — ~ \7> 7 \ > -D2 



2 K2(X 2 — X 1 )J ~ \2 K 2 (x 2 ~X 1 



2 



If we make the F-homotopic transformation 

H 3 (z) =z s i(z-lf*H 3 (z) 
with the requirement that of = —Bj then fH6|) reads 

+ $ s (z)#» + s (»tf s (z) = (47) 

with 

- l + 2ai 1 + 2S 2 ~ (ji + A 2 )z + « 2 + gi(l + 2a 2 ) - ji 

Since A 2 = -Ai it follows that (jUj) is an HYE. 

VI. THE METRIC : a = 1, c = >/2, 6 = k = I = , 7 = vr/4 

In this case we have 

t\ = 0, e 2 = -— , m(x) = —V2 x 2 , p(w) = 0, T(w, x) = wx + 1, 
3 3 

fW 2 (w) = fn^ n = h - Wi), X 2 (x) = g^ + Ax 3 - f 2 x 2 + g x x + g . 

n=0 i=l 

Taking into account that ([3]), and (jl]) imply in the present framework that 

9i = fa, 9o = 0, 
the function X 2 (x) can be written as follows 

3 

X 2 (x) = x (g^x 3 + fix 2 - f 2 x + f 3 ) = g A x JJ(a; - Xi). 

i=l 

Moreover, for i,j — 1, • • • , 3 let Wi and Xi denote the i-th root of the polynomial equations 
fW 2 {w) = 0, and X 2 (x) = 0, respectively. Throughout this section we shall assume that 
Wi 7^ Wj, and X{ 7^ Xj for every i ^ j. Furthermore, in the present case (jSJ) reads 

9* = - f ^ + fo 
19 



Finally, the roots of X 2 (x) = satisfy the following useful relations 

3 r 3 



x>=~' n 

i=l ^ 4 i=l 



A 

9 a' 



Making use of the expression for \l/ 2 given by (2.10g) in Kamran, and McLenaghan (1984), 
the functions f s (w), and g s (x) entering, respectively in (JSJ), and ([HD are computed to be 



f s (w) = -f 3 (l-\s\)(l-2\s\)w, 

r , w , „ (1- |s|)(l-2|s|) Tc/X 2 
2 v / 2re(s)(l - 2|s|)x - - — — 



dx x 



X 2 + f 3 



Hence, the equation for G s (w) becomes 

A ( f w 2 ^ 

dw \ dw 



with 



\ An 

j+T q ^ + Q s (w)G s = 

I-/ 2 



iV2q 



1 + / 2 



(48) 



Q s (w) := X s + / 3 (1 - \s\)(l - 2\s\)w + f\s\ (WW" + (1 - \s\)W' 2 ) - 



. r- W_ 2fq 2 1 

ivzsq w + (i + f 2 ) 2 w 2 - 

Equation fj48|) can be further simplified. To this aim let us transform G s (w) according to 

G s (w) = e h ^if s (w). 
If we require that h' = —T q /(2fW 2 ) we obtain the following equation for (p s (w) 



— I fW 2 ^ ] + R s (w) ^ = 



d_ 

dw 



dw 



(49) 



where 



R s (w) := X s + / 3 (1 - |s|)(l - 2\s\)w + f\s\ fww" + (1 - \s\)W'^ - 



r W' q 
iV2sq—— + — — 



W 2fW 2 ' 



For notational purposes let us define constants 

3 



Ci 1 ■= fsYliwi-wj), i = 1,2,3. 
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By means of the homographic substitution 



w — W\ 



(50) 



w 2 - Wi 

mapping the points wi, w 2 , W3, 00 to 0,1,25,00 with zs '■— (11*3 — wx)/(w 2 — w±), equation 
191) becomes 



d 2 y s . „, .dip, ~ , , 

— + P(z)— + Q s (z)<p s = 



(51) 



with 



p (2) = I + ^ + 1 



where for i = 1, 2, 3 



z z — 1 z — Zs 
B\ { B 2 B 3 



+ 



+ 



Ai A 2 A 3 
+ — + — — + 



z 2 (z — l) 2 (z — Zs) 2 z z — 1 z — Zs 



Bi = - ( - + ) , Ai = -(wx - w 2 )cig s (wi), 



Q s (Wi) = X s + f 3 

If we make the F-homotopic transformation 



C!-!H)(i.-2|H^-(^(2-| S |)-g 2 c 2 



/ 



Wj — 2w 



3=1 



p.(;0 = ^(z-lr(s-zsr&(s) 



and require that a 2 = —Bi, then (15T|) becomes 



Ci2 



with 



s l + 2ai l + 2o 2 l+2a 3 - , . A x A 2 A 3 

Z Z — 1 -2 — -25 ^ -2 — 1 -2 — Zs 



Finally, it can be checked that Ax + A 2 + A3 = with 

«i + 0:3(1 + 2ai) 



Ax = Ax — ax — a 2 (l + 2otx) 



zs 

a 2 + a 3 (l + 2a 2 ) 



A 2 = A 2 + ax + a 2 (l + 2a x ) 

z s — 1 

a 2 + a 3 (l + 2a 2 ) ax + a 3 (l + 2«i) 



v4 3 = A 3 + 



Z S -1 



+ 



^5 



(52) 
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Hence, (1521 reduces to an Heun equation with terms a(3 and accessory parameter qA given 
by (fT9l) . Concerning (Q we find the following equation for H s (x), namely 



^ ( X 2 ^ | + 



d 

dx 



dx 



-X s + AV2rsx + \s\ (XX" - (\s\ - 1)X' 2N ) - 2v / 2rsx 2 ^ r - 2r 2 ^-+ 
V / X X 2 

H s = 0. (53) 



X' 2 X 2 f 
(l-|s|)(l-2H) | 2— +' n 



ar x 



Let us introduce the following notation 



Cl 1 =5-4^1^2, C 2 1 =^4«l(jCl ~X 2 ), C 3 1 = g±X 2 {x 2 - Zl). 



By means of the homographic transformation 

X Xi — £3 

2 = 

x — X 3 Xi 

mapping the points 0, Xi, x 2 , x 3 , 00 to 0, 1, xs, 00, Xoo with 



X\ X3 X 2 Xqo 

, " — 



Xi 



x 2 - £3 



our equation (153]) becomes 



d 2 Hc „ , . <ii^ s . , . TT 

-— i + <p s ^ -r 2 + = 

a;<r dz 



with 



1 1 1 

+ 



Z Z — 1 -2 — Xg Z — Xqo 

#i B 2 B 3 2 



(54) 



(55) 



(56) 



Ai A 2 A 3 A r 
+ — + + — + 



Z 2 (z — l) 2 (z — Xs) 2 (z — Xqo) 2 z z — 1 Z — Xs z — x 



where 



A r 



Xi+ x 2 - 3x 2 

•£3X00 



B\ = , Bo 

4 



Xi — X3 
C2X1 ^jj ~ C3X2 ^jjj 



2 x 2 -x 3 . 



X3X5 



^1 = -— S > ^2 = — Qs , A 3 = — g s , g s = A s + — \s\A, 



0l I = X s + ^\s\B + 2V2rs^- 
2 xi - x 3 



1 - |s|)(l - 2|s|) [5f 4 (xi - x 3 )(xi - x 2 ) + / 3 xi] 



A 2 -~2 ^1^2 

4r c 2 # 4 - 



xi - x 3 



22 



= A s + %\s\C + 2V2rs^^- - (1 - - 2|s|) 
2 a; 2 - x 3 



g A (x 2 - x 3 )(x 2 - xi) - — 
4r 2 %g 4 %lA 



x 2 - x 3 ' 



A := (\s\ - 2)x 3 (xi + x 2 ) + 2(|s| - l)xix 2 , 

B := 2{2\s\ - 3)xj + (4 - 3\s\)x 1 x 2 + (M - 2)s 3 (x 2 - 2x0, 

C := 2(2|s| - 3)^ + (4 - 3|s|)xix 2 + (|s| - 2)x 3 (xi - 2x 2 ). 

If we make the F-homotopic transformation 

H s (z) = z Sl (z - lf 2 (z - x s f 3 (z - Xoo )H s (z) 

together with the requirement that af = —Bi, then (156|) becomes 

^ + y s (z)^ + & s (z)H s = (57) 

with 

fa, , 1 + 2ai 1 + 2S 2 1 + 2S 3 - 2li , 2t 2 2l 3 



z z — lz — xg z z — 1 Z — Xs 

Now, employing the expressions for the coefficients Ai with i — 1, 2, 3 it can be checked that 

• -, ■ , \Xqo Xoo Xs Xqq x 1 

i=l i=l ' 

Hence, (1571) becomes an Heun equation with 

a(3 = x s 2t 3 + 2l 2 , q A = -x s $li 

where again the 2lj's can be directly obtained from the formulae (J20l) -( j22l for the coefficients 
Ai by means of the formal substitutions Ai — > A* — > Ai, ai — > a«, zs — > xs and 

Zqo ^ Xqq . 

VII. THE METRIC <r 00 : b = £ = I, a = c = k = j = 
In this case we have 

e\ = 1 = p(w) = T(w, x), e 2 = = m(x), 

2 

fW\w) = ^2fnW n = f 2 (w- Wl )(w-W 2 ), 
n=0 

X 2 (x) = —Ax 2 + giX + g = —A(x — xi)(x — x 2 ). 
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Throughout this section we shall assume that w\ 7^ w 2 , and x\ 7^ x 2 . Making use of the 
expression for ^ 2 given by (2.11g) in Kamran, and McLenaghan (1984), the functions f s (w), 
and g s (x) entering, respectively in ([H]), and ^ are computed to be 

/» = |a(1-|s|)(1-2H), g s (x) = 0. 
The equation for G s {w) becomes 

^(^ 2 §) +P § + *WG^0 (58, 
with f := 2ir(f 2 - l)/(f 2 + 1), and 



2 W 
Q s (w) := X s - -A(l - |s|)(l - 2|s|) + 2irs— + 



f\s\ (WW" + (l-\s\)(W') 2 )+j^j- 2 



Ar 2 f 1 
2 



By means of the transformation G s {w) = e h ^ '(f s (w) with the requirement that ti 
—Y/(2fW 2 ) we obtain the following equation for (p s (w) 



^( /r S)^ sW ^ = ° (59) 



where 

^ 2 W' / 2\ r 2 

Q s (w) := A s - -A(l - |s|)(l - 2|s|) + 2trs— + f\s\ [WW" + (1 - \s\){W') J + 



W J 1 1 V v / fW r 

By means of the variable transformation z = (w — W\)/(w 2 — w\) mapping the points 
Wi, w 2 , 00 to 0, 1, 00 equation (1591) becomes 

^ + + Qs(z)<p. = (60) 

with 

P(z) = - + -, Q s {z) = — + - + — + 



z z — 1' Z 2 (2 — l) 2 z z — 1 

where 



2 f2(w 2 -w 1 )J ' \2 h{w 2 -W\) i 

.-H)(l-2| a |) | |a|(|s|-2) | 2r 

h 2 fi(w 2 - toi) : 



4 _ A,-(2/3)A(l-H)(l-2| 8 |) 1^1(1^1 — 2) 2r 2 _ 

A l _ ; 1 ^ r -pr, 777, /i 2 — —Al- 
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If we make the F-homotopic transformation 

and require that af = —Bi for i = 1, 2, then (1601) becomes 



with 



+ P s (z)— + QsUjy^ = 



dz 2 



dz 



(61) 



g f ^ _ 1 + 2a t | 1 + 2a 2 ^ _ A 2 )z + a 2 + a^l + 2a 2 ) - A 

Z Z — 1 Z(2 — 1) 

Since A 2 = —A\ it follows that (IBTj) is an HYE. Concerning ([9]) we find the following equation 
for H s (x), namely 

d ( „ndff«\ [ . , , / // , „/2\ _ X g 2 



H„ = 0. 



(62) 



dx \ dx J [ s ' 1 1 \ ' U ^ ^ X X 2 

By means of the transformation z = (x — Xi)(x 2 — Xi) mapping the points x±, x 2 , oo to 0, 1, oo 
our equation (|62|) becomes 



with 



where 



d 2 H„ _ . . dH* , . 

—l.+y a ( z )—!L+£l a ( z )H e = 
dz z dz 



™ , ^ 1 1 ^ , n B\ B 2 A x A 2 

y a (z) = - + -, Q s (z) = -± + ~ — + — + ' 



(63) 



z z — 1 



z 2 (z — l) 2 Z 2 — 1 



#1 
At 



s 

■x + 



Q 



2 A(x 2 -x0; ' 
A. |s|(|s|-2) 2g 2 
A " 



So = - - 



9 



2 A 2 (x 2 -x0 2 ' 
Finally, with the F-homotopic transformation 

H s (z) = z s ^z-irH s (z), 

and the requirement that of = — for z = 1, 2 (|63|) becomes 

d 2 H * -r- . . (ii/o t , . ~ 
_i +!p , W _£ +iW = 



2 A(x 2 -xi) 
1 2 = -A- 



(64) 



with 



^ _ 1 + 257! | 1 + 2a 2 ^ {A 1 ^A 2 )_z + a 2 + a 1 {l + 2a 2 )-A 1 



z z — 1 ' z(z — 1) 

Taking into account that A 2 = —A\ it results that fl64|) is an HYE. 
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